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Turbulent viscosity in clumpy accretion disks II 
supernova driven turbulence in the Galaxy 
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Abstract. An analytical model for a turbulent clumpy gas disk is presented where turbulence is maintained by the 
energy input due to supernovae. Expressions for the disk parameters, global filling factors, molecular fractions, and 
star formation rates are given as functions of the Toomre parameter Q, the ratio between the cloud size and the 
turbulent driving length scale S, the mass accretion rate within the disk M, the constant of molecule formation a, 
the disk radius, the angular velocity, and its radial derivative. Two different cases are investigated: a dominating 
stellar disk and a self-gravitating gas disk in z direction. The turbulent driving wavelength is determined in a first 
approach by energy flux conservation, i.e. the supernovae energy input is transported by turbulence to smaller 
scales where it is dissipated. The results are compared to those of a fully gravitational model. For Q = 1 and 8 — 1 
both models are consistent with each other. In a second approach the driving length scale is directly determined 
by the size of supernovae remnants. Both models are applied to the Galaxy and can reproduce its integrated and 
local gas properties. The influence of thermal and magnetic pressure on the disk structure is investigated. We 
infer Q ~ 1 and M ~ 0.05 - 0.1 M Q yr _1 for the Galaxy. 

Key words. ISM: clouds - ISM: structure - Galaxy: structure - Galaxies: ISM 



1. Introduction 

The interstellar medium (ISM) is of turbulent nature. Its 
structure is usually described as hierarchical (Scalo 1985) 
over length scales of several magnitudes up to ~100 pc. 
The neutral phase of the ISM is not uniform but of fractal 
. ' na t ure (Elmegreen & Falgarone 1996). 

In a widely accepted picture (see, e.g. Kulkarni & 
Q ■ Heiles 1988, Spitzer 1990, McKee 1995), the ISM consists 
5^ , of 5 different phases that are listed in Table ^ About 
80% of the total gas mass is neutral and 50% is in form of 
clouds or filaments. The scale height of the cold material is 
about 4 times smaller than that of the warm gas. The hot 
medium is heated and ionized by direct thermal energy in- 
put due to supernova (SN) explosions. The warm medium 
is heated and ionized by SN remnants, stellar winds and 
radiation. The main energy sources driving interstellar 
turbulence and thus causing the multiphase structure of 
the ISM are (i) SN explosions (see, e.g., Ruzmaikin et al. 
1988) or (ii) the galaxy's gravitational potential and local 
gravitational instabilities (Wada et al. 2002). 

Several attempts have been made to include SN explo- 
sions into hydrodynamical simulations: 

Rosen & Bregman (1995) computed two-dimensional 
simulations with two cospatial fluids representing stars 
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and gas. Star formation, stellar mass loss, stellar winds 
and SN were included. These simulations created a three- 
phase medium with filaments of dense, cold and warm gas 
surrounding bubbles of hot gas. 

Korpi et al. (1999) simulated the local ISM heated and 
stirred by SNe using a three-dimensional, non-ideal MHD 
model in a box of ~1 kpc. They used a prescription for 
the SN rate deduced from observations. These simulations 
showed a stationary multicomponent structure of the gas 
in a state of developed turbulence. They found a turbulent 
cell size of 60 pc in the warm phase. 

Gazol-Patino & Passot (1999) studied the role of su- 
perbubbles in the evolution of the ISM at the kilopar- 
sec scale. The model incorporated the fully compress- 
ible MHD equations, including parameterized cooling and 
heating, the Coriolis force, shear, self-gravitation, and en- 
ergy input by SNe in 2.5 dimensions. The SN rate was 
determined by the gas flow itself. Their simulations also 
lead to a state of highly compressible turbulence. The in- 
clusion of SNe into the model gave rise to superbubbles 
whose presence leads to a more intermittent cycle of star 
formation. 

Wada & Norman (2001) made high-resolution two- 
dimensional hydrodynamic simulations of the ISM includ- 
ing star formation and energy input by SNe. In their sim- 
ulations a globally stable multiphase ISM formed. They 
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T(K) 


n (cm 3 ) 


< n > (cm -3 ) 


?; rms (kms 


$v 


H (pc) 


M/Mtot (%) 


hot ionized 


~ 10 b 




0.002 




0.5 


3000 


4 


warm ionized 


~ 8000 


0.3-10 


0.025 


~10 


0.2 


900 


14 


warm neutral 


- 8000 


0.1-10 


0.1 


~10 


0.3 


400 


31 


cold neutral 


~ 100 


10-1000 


0.3 


~6 


0.02 


140 


25 


molecular 


~ 10 


> 100 


0.6 


~6 


0.001 


70 


26 



Table 1. The properties of the different phases of the ISM (T: temperature, n: density, < n >:mean density, v. 
dispersion velocity, $v : volume filling factor, H: scale height, M: percentage of the total gas mass (from Boulares & 
Cox 1990). 



found an energy spectrum of E{k) oc k 2 and a turbulent 
length scale of 20-100 pc. 

The influence of SN explosions on the ISM was studied 
by de Avillez & Mac Low (2001) with a three-dimensional 
hydrodynamical code using adaptive mesh refinement. SN 
explosions were set up at a rate comparable to observa- 
tions. They showed that hot gas that is not evacuated 
through chimneys expands into the cooler gas of the thick 
disk building mushroom-shaped structures. 

Even though it is clear that SN explosions represent 
the dominant energy input in galaxies with a normal star 
formation rate, all these simulations show how SNe can 
effectively drive ISM turbulence in detail. 

In order to form massive stars which have a SN ex- 
plosion at the end of their lifetime, self-gravitation is the 
starting point. Self-gravitation drives star formation and 
subsequently SN explosions. 

In a previous paper (Vollmer & Beckert 2002) we have 
elaborated a model of an equilibrium disk where turbu- 
lence is generated by instabilities involving self-gravitation 
and maintained by the energy input from differential ro- 
tation and mass transfer. The model only involves gravity. 
Its application to the Galaxy showed good agreement with 
observations. 

In this article we extend our previous model by includ- 
ing SNe as an additional energy input to drive ISM turbu- 
lence. We give the basic picture in Sect. [2] The equations 
are presented in Sect. [21 Two different models with differ- 
ent vertical pressure equilibria are analyzed: (i) a model 
of a self-gravitating gas disk in z direction (SGZ model) 
and (ii) a model with a dominating stellar disk mass (DSD 
model). We introduce an energy flux conservation equa- 
tion that relates the energy flux transported by turbulence 
to the energy input due to SN explosions. In Sect.0]we give 
analytical expressions for the disk properties and show re- 
sults for characteristic values of the input parameters. The 
SGZ model is compared in Sect.[5]to the fully gravitational 
model of Vollmer & Beckert (2002). In Sect.Elthe energy 
flux conservation is replaced by an equation relating the 
turbulent driving length scale directly to the size of a SN 
remnant (SNR model). We apply the SGZ and SNR model 
to the Galaxy in Sect.[7| In Sect.|S|a realistic gravitational 
potential for the galaxy is introduced and the inclusion 
of magnetic and thermal pressure into the models is dis- 
cussed. Sect. ITfil investigates disks with Q > 1 and Q < 1. 
The summary and conclusions are given in Sect. 1111 



2. The basic picture 

We consider the warm, cold, and molecular phases of the 
ISM as one gas that can undergo phase changes accord- 
ing to internal and external conditions. Internal conditions 
are pressure, metalicity, molecule formation, and turbu- 
lent time scales; external conditions are stellar radiation 
field (UV, X-rays, cosmic rays), stellar winds, and SN 
explosions with subsequent shock formation. In this ap- 
proach we neglect the thermal balance due to radiative 
heating and cooling of the ISM. 

The ISM is assumed to be turbulent. Turbulence is 
driven by SNe energy input at a characteristic length scale 
^driv This length scale might be identified as the charac- 
teristic length scale of a SN bubble or alternatively that 
of the interaction of SN bubbles. The disk scale height 
is determined unambiguously by the turbulent pressure 
Pturb = P w ?urb' wnere P is the average density and Vturb 
the turbulent velocity in the disk. We thus neglect ther- 
mal, cosmic ray, and magnetic pressure. The interaction 
of SN bubbles leads to viscous transport of angular mo- 
mentum. We assume that the energy input due to SNe is 
dissipated by viscous heating. 

Thus, we consider a gaseous turbulent accretion disk 
in a given gravitational potential $ which gives rise to an 

angular velocity ft = J R^ 1 ^- The Toomre parameter 
(Toomre 1964) is treated as a constant but free parameter 
for the whole disk 

Q = VJ ^, (i) 

with the restriction Q > 1, where G is the gravitational 
constant, E the average gas surface density, and k the local 
epicyclic frequency. 

3. The equations 

3.1. The volume filling factor 

We compare the crossing time of a turbulent cloud to the 
gravitational free fall time in order to derive an expres- 
sion for the volume filling factor (j>y. The characteristic 
turbulent time scale of clouds whose size is a factor 6~ !+ D 

_ 3 

smaller than the driving length scale l c \ = 5 1 + D Zdriv is 

tl = S'i'^ Zdriv/Vturb , (2) 
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where D is the fractal dimension (see, e.g., Frisch 1995). 
The local gravitational free fall time is given by 



3tt 



32Gp c 



(3) 



where p c \ is the density of a single cloud, which is related to 
the overall disk density p by the volume filling factor cf>y: 
Pel = <ihj X P- The clouds become self- gravitating for t\ = tg. 
We assume D = 2 for a compressible, selfgravitating fluid, 
which is close to the findings of Elmegreen & Falgarone 
(1996). 



3.2. The viscosity prescription 

In general the viscosity v is defined as the product of 
the characteristic velocity v and the characteristic length 
scale / of the system: v = vl. We assume that turbulence 
becomes intermittent due to self-gravitation of cold gas 
clouds. There are two simple frameworks for understand- 
ing intermittency: (i) based on velocity (ii) based on dis- 
sipation. 

Within the framework of intermittency based on dis- 
sipation (Frisch 1995), the local spatial average of the en- 
ergy dissipation rate per mass unit e, which is constant 
over all length scales in the case of a Kolmogorov-like tur- 
bulence, depends explicitly on the length scale 

3-D 

(4) 



"^turbAdr 



where D is the fractal dimension (D=3 for a Kolmogorov- 
like turbulence), l drlv is the driving length scale. This can 
be achieved if there is an energy sink in the turbulent cas- 
cade. We identify the self-gravitation of cold gas clouds as 
this energy sink, i.e. during the contraction of a cloud kine- 
matic energy is transformed into heat and radiated away 
so that the cloud can further shrink. Thus self-gravity 
gives rise to intermittency in the ISM turbulence. 

Within the framework of intermittency based on ve- 
locity, intermittency can be understood as a decreasing 
volume filling factor of turbulent eddies with decreasing 
length scale. This can also be achieved by self-gravitation, 
because self-gravitating clouds are smaller than diffuse 
clouds of the same mass. 

In the case of intermittent turbulence the turbulent 
viscosity v — fturt^driv overestimates the true viscosity 
and has to be decreased by a factor 7: 



V = 7 w turb^dr 



(5) 



This factor is related to self-gravitation, i.e. the free fall 
time of the disk in z direction. Based on these arguments 
we make the hypothesis that the characteristic velocity is 
determined by the disk height H and the averaged local 
free fall time t^: v = H/t^. The appropriate length scale 
is the turbulent driving length scale /driv 
Thus we obtain: 



V — ^driv ■ 
% 



(6) 



The free fall time is given by = \J (3n)/(32Gp), where p 
is the averaged midplane density of the disk. Using the lo- 
cal free fall time (Eq. ©) with p c \ = (p^p, where 4>y is the 
volume filling factor of self-gravitating clouds fSect. l3~Tj) . 
The characteristic velocity can then be written as 

VWWturb • 



Inserting this velocity into Eq. JHJ) gives 

V = HyQ>y~6v turh . 

Thus, the factor 7 that is due to intermittency is 
v<5 ■ 



H 



7 



dri' 



(7) 



(8) 



(9) 



We will call it the viscosity intermittence factor. 
3.3. Angular momentum equation 

The viscosity prescription given above is essential for mass 
accretion and the transport of angular momentum. In a 
steady state accretion disk, the mass accretion rate is 

M = 2nRZ(-v Iad ) , (10) 

where w ra( j is the radial velocity. The angular momentum 
equation can be integrated giving 



2ttR K dR J 
Furthermore, we use 



(11) 



Y> = pH (12) 

for the surface density of the disk, where p is the average 
mass density at the midplane. 

3.4. The vertical pressure equilibrium 

We assume that the only pressure which counterbalances 
gravitation in the vertical direction is the turbulent pres- 
sure pturb = P^turb ' We distinguish two cases for the grav- 
itational force density in the vertical z direction: 

1. p < p* and Md(R) < M*(R), where p*/M*(R) is 
the stellar central density/disk mass within a radius R 
(dominating stellar disk mass); 

2. p > p* and 0.5 (H/R) M{R) < M d (R) < M(R), where 
M(R) is the total mass enclosed within a radius R 
(self-gravitating gas disk in z direction). 

In the following we call the model of a dominating stellar 
disk mass DSD model and the model of a self-gravitating 
gas disk in z direction SGZ model. The hydrostatic equilib- 
rium in the vertical direction implies that the gravitational 
force is balanced by the turbulent pressure. 

For the two cases the hydrostatic equilibrium has the 
following forms: 

1. pturb = £9$/<9z , (13) 
where $ is the gravitational potential of the disk; 

2. Pturb = ttGS 2 (14) 
(Paczynski 1978). 
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3.5. Global gravitational stability in z direction 

The basic principles underlying the gravitational instabil- 
ity of a thin rotating disk can be found in Toomre (1964). 
A gaseous disk is locally stable to axisymmetric perturba- 
tions, if 



7T(j- Zj 



(15) 



where n = y R^j^ + 4f2 2 is the epicylcic frequency. Since 
in general Q < n < 2f2, we will use the following equation: 

^ V t ur b H 



vrGS 



(16) 



Multiplying the numerator and the denominator of the 
right hand side by R 2 gives 



^turb M tot 



M. 



(17) 



gas 



where M< 



tot /gas 



is the total enclosed mass and the total 



enclosed gas mass at radius R. Thus for a given velocity 
dispersion Q _1 is proportional to the ratio of gas mass to 
total mass. 



3tt/(32G( 



Pel 



500 enr 



r )) ~ 2 Myr . 



(23) 



Since <tmb = Sig, this leads to <5=5, which is precisely 
the value we find in order to fit the gas properties of the 
Galaxy (Sect. CJ). 

3.7. Energy flux conservation 

Here we assume that the energy input into the ISM due 
to SNe is transported practically without loss from the 
driving length scale to the dissipative length scale where 
it is radiated away. Mass is stored into self-gravitating 
clouds and eventually used for star formation. The energy 
per unit time which is transferred by turbulence is 



E 



P v \ "turbAdriv d ^ . 



(24) 



where the integration is taken over the volume J dV = 
V = AH (Landau & Lifschitz, 1959). Thus, the energy 
flux per unit time and unit area is 



AE 
AtAA 



= -Si/ 



"turb 
I 2 



(25) 



3.6. The star formation rate 

It is generally accepted that the star formation rate is pro- 
portional to the mean density of the disk and the inverse 
of the characteristic time scale for the cloud collapse, i.e. 
the non-averaged local free fall time t s (see Sect. 13. l|) : 

P 

P* oc -j- . (18) 



Since tg cx p 2 (see Eq. (01) this corresponds to a Schmidt 



law of the form oc p? . The factor of proportionality is 
given by the probability to find a self-gravitating cloud, 
i.e. the volume filling factor (f>y. Thus, the star formation 
rate is given by 



P_ 

'ir 



(19) 



Furthermore, we assume that stars are only born in the 
midplanc of the disk in regions that have the size of the 
turbulent driving length scale Zdriv, because the clouds 
can collapse only within the turbulent time scale tturb = 



^driv/^turb- We thus obtain 



driv 



(20) 



for the mass surface density turned into stars. The vol- 
ume filling factor is defined such that t\ = t l s = <5 -1 £turb 
(Sect. EH). Combining Eq. (UHJ and Eq. (JSOJ we get 



Hurb 

A simple estimate gives: 

^driv \ / ^turb 



^turb 



100 pc y V 10 kms 



-) 1 - 10 Myr and 



(21) 



(22) 



This energy loss is balanced by the SNe energy input 



SN 



E = 



(26) 



We connect this energy input directly to the star forma- 
tion rate (Eq. (|20Jl ). With the assumption of a constant 
initial mass function (IMF) independent of the environ- 
ment one can write 



E$n 
~AA 



(27) 



The factor of proportionality £ relates the local SN en- 
ergy input to the local star formation rate and is taken to 
be independent of the radius. It can be normalized with 
Galactic observations by integrating over the Galactic 
disk: The number of SN per time is proportional to the 
Galactic integrated star formation rate: 

iVsN = f M» . (28) 

The energy input per surface area and time in Eq. I|25|) is 



AE S j?iVgN ^SN^-^* 



AAAt 



J SN J 

AA 



AA 



(29) 



where E^ is the kinetic energy input from a single SN. 
Thus the energy released into the ISM per mass turned 
into stars is 



C NsN pkin 



(30) 



Thornton et al. (1998) have shown by modeling SN 
explosions in different environments that the kinetic en- 
ergy of the remnants is ~10% of the total SN energy E<^ 
irrespective of the density and metalicity of the ambi- 
ent medium. The SN energy input into the ISM is thus 



B. Vollmer and T. Beckert: Turbulent viscosity in clumpy accretion disks II 



5 



pkin 



10 ergs. The integrated number of SNe type II 
in the Galaxy is taken to be Asn ~ 1/40 yr _1 in ac- 
cordance with Tammann et al. (1994). The Galactic star 
formation rate is taken to be M*=3 M Q yr -1 (Prantzos & 
Aubert 1995). This leads to f ~ 4.6 10~ 8 (pc/yr) 2 . 



3.8. The molecular fraction 

In order to derive an expression for the molecular fraction 
of gas in the disk, we compare the crossing time of the 
turbulent layer iturb and the H-H2 transition time scale 
tn 2 = a$ v /p (Hollenbach & Tielens 1997). We define the 
molecular fraction here as / mo i = iturb / tn 2 . This allows us 
to calculate the molecular surface density and total mass 
that can be compared with observations. 



4. Results 

With £ fixed and with a given rotation curve our model has 
only three free parameters: the constant of the molecule 
formation rate a, the Toomre parameter Q and the mass 
accretion rate within the disk M. All disk properties can 
be expressed as functions of a, Q, 8, M, £ f2, dfl/dR, and 
R. We have solved the set of equations Eq. JHJ, Eq. (|1 1|) . 
Eq. H2J), Eq. JE)/©, Eq. (JTSJ, and Eq. J2HJ), for the 
two cases: 

1. dominating stellar disk (DSD) (Q > 1, SO and gas 
deficient galaxies), 

2. self-gravitating gas disk in z direction (SGZ) (Q ~ 1, 
spiral galaxies). 

We will use f2' = dVt/dR with il 1 < in the region of 
interest. 



4.1. Dominating stellar disk (DSD) 

In this case the vertical pressure equilibrium is given by 
Pgrav — ^d^/dz. We make the approximation d$/dz ~ 
Xw 2 ot ^ _1 leading to p grav = x^GS^E (see e.g. Binney & 
Tremaine 1987). The factor x allows a smooth transition 
between the SGZ model (Sect. l4~2|) and the DSD model. 

As in Sect. 14. 2l we first give the expression for H, idriv, 
p, and X as functions of fturb- 



Vrot 



(31) 



*driv = J\\x h GMQi v J b & R -i , (32) 



Wturb 



(33) 



S = 7r- 1 G- 1 Q- 1 t; tU rb^ . (34) 
The expression for the turbulent velocity dispersion is 
v turh = 0.87x™G*M*Q™6-™£™R™ilTs(-n')-* .(35) 



Inserting Eq. (JSSJl into Eq. and 01 leads to 

H = 0.75 x~*G*M*<2* ^R~^n~^(-Q')~^ ,(36) 

Zdriv = omx~*GHi?Qh^-^R~ i n-i(-n , )-' ! , (37) 

£ = 0.28x™G-§MiQ-^5-&(TER^QTi(-n')-i . (38) 



The viscosity u, the volume filling factor 0y, the star for- 
mation rate £*, and the molecular fraction / mo i then write 



v = o.58x _ *G*M*gT*jT5£-T5i?-Tf n-*(-n')-% , (39) 

<fa = 0.79x~*G*M*Q-ft(TTi£-ftirAftA(_r2'r* ,(40) 
£* = 0.25x#G-*M&Q-S$-A£-*i£&fi8(-n')-§ ,(41) 
J mo i = 0.37x^G-T§A/-^g-^£*i?T*fiTl(„^')^ Q ,-i .(42) 



Setting f2' ~ —tt/R and u ro t = fii? leads to the follow- 
ing simple radial dependences: Wturb = (x w rot) TC j H oc 

(XWrot)^"^" 1 , Zdriv OC (\ «rot ) ~ 5 f^ 1 , £ CX (x^rot) 15 ^, 
^ (X (x«rot)~^^\ 4>Y = (X V Tot)~™, £* OC (X^rot)^fi 2 , 

and / mol oc (xWrot) 36 "^- 

We found that for Q=l and x ~ 0-05 there is a 
smooth transition between the SGZ and DSD model, thus 
d&/dz ~ 0.05u 2 ot i? _1 . Fig. 2] shows the radial depen- 
dences of the above parameters for (i) a constant rota- 
tion curve and (ii) for a rising rotation curve. We have 
included an arbitrary rising rotation curve to show the 
general effects of (— O') on the disk properties. We chose 6 
in a way to fit the observed volume filling factor of molec- 
ular clouds (Table [TJ. For this model x = 0.05, Q = 10, 
6 = 5, M = 10- 1 Mgyr" 1 , and a = 3 10 7 yrMopc" 3 . 
The rising rotation curve leads to a bigger driving length 
scale, disk height, viscosity, and volume filling factor. On 
the other hand, it leads to a smaller surface density, den- 
sity, molecular fraction, and star formation rate in the in- 
ner disk than for a constant rotation curve. The total gas 
mass is M gas - 1.3/1.0 10 9 M Q , the total molecular mass 
is M gas ~ 1.7/0.5 10 8 Mq, and the total star formation 
is M* ~ 0.2/0.1 M Q yr _1 for the model with a flat and a 
rising rotation curve, respectively. This model applies to 
gas poor spiral and SO galaxies. 

For this model we have used £* = p» idriv This is only 
valid if H /idriv > 1 with 



H 



= 1.09x~ TO G"ToM-ToQ-io ( rTo ( f 



(43) 



For w rot =220 kms -1 , M = 10~ 2 Moyr" 1 , and x = 0.05 
this translates into the requirement Q < 1440. 

In the case of H/l^ lv < 1, we set £* = p* H. Whereas 
Eq. H31f) ~ Eq. (|34|l do not change, the turbulent velocity 
becomes 



v tmh = 0.87 x^G^M^Q*5~^£,^R^fl^ (-0' v ~~ 



.(44) 



The exponents are only slightly different from those of 
Eq. ESJ- 
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Fig. 1. Disk parameters for the case of the DSD model. 
Solid lines: constant rotation curve. Dashed line: rising 
rotation curve v YOt oc y/~R. 

4.2. Self-gravitating gas disk in z direction (SGZ) 

We first give the expressions for the disk height H : the 
turbulent driving length scale ^riv 5 the mean density in 
the disk plane p, and the disk surface density £ as func- 
tions of the Toomre parameter Q, the mass accretion rate 
M, the angular velocity f], the disk radius i?, the radial 
derivative of the angular velocity £7', and the turbulent 
velocity dispersion vturb- 



H 



/dri 



Qv tv 



n 



3 ttGMQ 
28Rv* wh (-Sl>) 



n 2 



ttGQ 2 

^turb 



(45) 
(46) 

(47) 
(48) 



nGQ 

The expression for the turbulent velocity dispersion is 
uturb = 0.82G&M&6-&£&R-&nfr(-Q')-fi: . (49) 



Since Q = (v tmh /v TOt )(M tot /M sas ) (Eq. JHJ) and the 
turbulent velocity dispersion Wturb is independent of Q, the 
Toomre parameter is only a measure of the ratio between 
the total enclosed and the total gas mass of the galaxy. 

Using Eq. © together with the pressure equilibrium in 
z direction the viscosity intermittence factor 7 of Eq. JSJ 
can be calculated 



7 



H 

^driv 



32 
3^ 



1 . 



(50) 



Thus, in the SGZ model the viscosity reads 

V = V^V Wturb H = Uturb 'driv ■ (51) 

Inserting Eq. gSJ into Eq. gSJ, gHJ), and (gHl leads 

to 

H = 0.82G&M*Q<r*f*ir*fr*(-n')~* , (52) 

ldriv = 0.71G*ilf*Q*&£-*Jr*n-*(-n / )~* , (53) 
S = 0.26G-ftM*Q- 1 5-n^Ji-AnTf(-fi')-ft . (54) 

For the viscosity v, the volume filling factor (j)y, the star 
formation rate £», and the molecular fraction / mo i one 
obtains 

v = o.6iG&MAQ<?Af _ A.R _ Afi-A(-fi / )-& , (55) 
V = o.8ig*m ^r^£-Air^rm(-iy)-TT , (56) 

£* = 0.21 G-nMirQ-^-nf fiJ-^off > ( 57 ) 

/ mol = 0.34G-TT M-TiQ-^TTR^nrTi-n'^a- 1 . (58) 

Setting fl' ~ —Sl/R leads to the following simple radial 
scalings: vtmb = const., H oc f2 _1 , idriv oc Q , E oc f2, 
oc f2 , (/>v = const., E* oc fi 2 , and / mo i oc f2. 

Fig. shows the radial dependences of these parame- 
ters for (i) a constant rotation curve and (ii) for a rising 
rotation curve. We have again included a rising rotation 
curve to show the general effects of (— f2') on the disk prop- 
erties. For this model Q = 1, 5 = 5, M = 10 _1 M yr _1 , 
and a = 3 10 7 yrM Q pc -3 in order to fit Galactic observa- 
tions. The rising rotation curve leads to a bigger driving 
length scale, disk height, viscosity, and volume filling fac- 
tor. On the other hand, it leads to a smaller surface den- 
sity, density molecular fraction, and star formation rate 
mainly in the inner disk. The fraction between the gas den- 
sity, and star formation rate for a constant and those for 
a rising rotation curve is about 5-10 in the central part of 
the galaxy. This fraction is about 3-5 for the surface den- 
sity, driving length scale, disk height, and viscosity. The 
difference in the velocity dispersion and the volume filling 
factor is 20% and 30% respectively. The total gas mass for 
both models is M gas ~ 6.8 10 9 M Q , the total molecular 
mass is M mo i ~ 3.0 10 9 M Q , and the total star formation 
is M* ~ 2.6 M©yr _1 in good agreement with Galactic ob- 
servations (for the gas see, e.g. Kulkarni & Heiles 1988; for 
the star formation rate see, e.g. Prantzos & Aubert 1995). 
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Fig. 2. Disk parameters for the SGZ model. Solid lines: 
disk with constant rotation curve. Dashed line: disk with 
rising rotation curve u ro t oc ^/R. 



5. Comparison with the fully gravitational model 
of Vollmer & Beckert (2002) 

In Vollmer & Beckert (2002) we have analytically investi- 
gated the equilibrium state of a turbulent clumpy gas disk 
that consists of distinct self-gravitating clouds, which are 
embedded in a low density medium, and evolve in the fixed 
gravitational potential of the galaxy. Gravitational cloud- 
cloud interactions in the disk give rise to an effective vis- 
cosity and allows the transport of angular momentum and 
mass in the gas disk. In this scenario turbulence is assumed 
to be generated by instabilities involving self-gravitation 
and to be maintained by the energy input, which is pro- 
vided by differential rotation of the disk and mass transfer 
to smaller galactic radii via cloud-cloud interactions. Both 
the energy source and dissipation process for turbulence 
in these disks are due to gravity. Therefore, we call this 
model the fully gravitational model. 



5.1. The viscosity prescription 

In the fully gravitational model the viscosity prescription 
is 



v = Re Vt 



urb 



H 



(59) 



In Vollmer & Beckert (2002) we have used a slightly dif- 
ferent expression for the energy flux transferred by tur- 
bulence to smaller length scales: E = p^v^ nl . h /ldn v = 
£w t 2 m . b /(in driv ) instead of E = S^ t 2 urb /^ r iv Thus we 
had already assumed /driv = H. This lead to exponents 
of Q that were twice that of the model using Eq. (|25|) . 
Therefore, the volume filling factor in the case of a tur- 
bulent non-Kolmogorov energy spectrum (D = 2) of the 
form E(k) oc fc~ 2 in the framework of Vollmer & Beckert 
(2002) using Eq. (gj is 



Q~ 2 Re 



(60) 



For Q = 1 one gets 4>y — Re 2 . In the present model 
including SNe, Eq. JHJ reads for the viscosity 



v = VrV fturb 5 H . 



(61) 



Thus, for Q = 1 and 6=1 the viscosity prescriptions of 
both models are equivalent. 

5.2. The star formation rate 

In Vollmer & Beckert (2002) we suggested a star formation 
law of the form 



Re-^fl 



(62) 



This has to be compared with the star formation rate in 
the present model, which reads (see Sect. 13. 6|) . 



— V^V P^driv (iff ) 1 



(63) 



In the case of the SGZ model p = (tt G Q 2 )" 1 ^ (Eq. gZJ). 
Thus, for Q = 1, (tg) _1 — fi leading to a star formation 
rate of the form 



(64) 



Since in the fully gravitational model Idi-w = H, the pre- 
scription for the star formation rate of both models are 
equivalent for Q = 1 disks. 



5.3. Energy flux conservation 

The main difference between the fully gravitational model 
and the model including SNe lies in the energy flux con- 
servation equation. The energy transferred by turbulence 
to smaller length scales is balanced by the energy input 
through SNe in our present model, whereas it is balanced 
by the differential rotation and mass inflow in the fully 
gravitational model (Vollmer & Beckert 2002). 
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5.4. Results 

If one approximates 0! ~ —fl/R, the radial scalings of 
all disk properties is the same for both SGZ models, i.e. 
fturb = Const., H OC fi~ , Zdriv oc fl^ 1 , E oc f2, ^ oc f2 _ , 
0v = const., E* oc f2 2 , and / mo i oc fi. 

We therefore conclude that the fully gravitational 
model (Vollmer & Beckert 2002) and the here presented 
model including SNe lead to the same radial dependences 
of the disk properties in the case of Q = 1 with Re = (f>y . 
A quantitative comparison will be made in Sect.0 

For the comparison with previous models of turbulent, 
self-gravitating gas disks we refer to Sect. 5.1 in Vollmer 
& Beckert (2002). 

6. SN remnant size as driving wavelength (SNR) 

In Sect. 13.71 the energy flux conservation equation 
determines the turbulent driving wavelength ^driv 
Alternatively, one can assume that only the interaction of 
SN bubbles leads to an effective turbulent viscosity that 
transports angular momentum. In this case the driving 
wavelength is twice the radius of the SN remnant at a 
time when the shock thermalizes. In the following we will 
call this model the SNR model. Following Dorfi (1993) 
the final radius of a SN remnant is given by 



Rsnr — IsnE ; 



51 "o 



pc 



(65) 



where ^sn is the characteristic radius for no=l cm 3 , 
E51 — Ej (10 51 erg), and no is the average local density in 
cm- 3 . We set ; S n=64 pc (Dorfi 1993) and replace Eq. 
by IdHv = 2i?sNR- This results in 

wturb = 0.56 Got Q*£- 



IsnMzR 



H = 0.56GOTQ7o£ 51 ™; sN 2 M 2^- 5 ft-35(_ry)-2 

13 26 . ■ 1- . 26 

P ~ nGQ 2 ' 

24 61 — • 1 1 48 

E = 0.18 G-^Q-^E'^l^M iR-in^(-n')- 

24 61 l_L . 1 i i ^ 1 

i/ = 0.89GMQ7o£'™/| N M2 J R-2n-3i(-fi')-2 , 
$v = 8.05 G^Q-^S^E^l^M^RCl-^i-n') 



(66) 
(67) 
(68) 

(69) 

(70) 
(71) 
(72) 



E« = 1.44G~TQ-fitf-i^Mi| N J|i-l J Rin¥(_n / )* , (73) 



/ mo i = 0.11 G~-s Q~ ™5E 51 70 l s ^M 2 R-2Q-s{-n')- 2 a' 1 

The intermittence factor 7 (Eq. ©) is in this case 7 ~ 1 
as for the SGZ model. 

Setting fl' ~ —Cl/R leads to the following radial scal- 
ings: v tuIh oc rr 9 / 70 , h oc rr 79 / 70 , z driv cx rr 26 / 35 
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Fig. 3. Disk parameters for the case of SNR model. Solid 
lines: disk with constant rotation curve. Dashed line: disk 
with rising rotation curve v ro t oc yf~R. 



and / mo i oc ft 43 / 70 . The exponents are quite close to those 
of the SGZ model. 

We assume M = 10 _1 Moyr -1 , Q = 1, S = 5, and 
a = 2 10 7 yrM©pc -3 . Fig. shows these parameters for 
a constant and a rising rotation curve. In comparison to 
the SGZ model, the turbulent velocity dispersion and the 
volume filling factor are not constant with radius. For a 
constant rotation curve the latter increases by a factor 4 
from the outer disk to the center of the galaxy. All other 
quantities are similar to those of the SGZ model. 

The total gas mass for both rotation curves is M gas ~ 
6.3 10 9 M and the total molecular gas mass is M mo i ~ 
j-2^10 8 Mq. The total star formation rate is M* ~ 
0.8 M Q yr _1 for the rising rotation curve and M* ~ 
5.7 M0yr _1 for the constant rotation curve. 

The star formation rate for a constant rotation curve 
is a factor 2 higher than the observed Galactic value. 
Alternatively, if one drops Eq. (|20(l in the initial equa- 
tions, one can deduce the star formation rate a posteriori 
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with the help of the energy flux conservation equation 
(Eq. (El): 

S, = r 1 S^(^) 2 (75) 

tdriv 

This is shown for the constant rotation curve as a dot- 
ted line in the panel of the star formation rate in Fig. [3] 
Since Y,v ~ M/(3ir) and v tln h = const. Eq. I|75|) yields 
£* oc l^ iv . Now, with ^driv oc p" 13 / 35 (Eq. p cx SI 2 

(Eq. dHHJ), and E cx SI 61 / 70 one obtains cx p ' 75 cx 
Vt 1 - 5 cx S 1 ' 7 , which is the observed Schmidt law. In this 
case the characteristic timescale for star formation is no 
longer the free fall time but it is proportional to the vis- 
cous timescale. The total star formation rate in this case 
is M* ~ 1.8 M yr -1 . 

A possible remedy would be to use the star forma- 
tion prescription of Eq. i|75|) an d assume a smaller £, 
i.e. a smaller fraction of the total SN energy that goes 
into turbulent motions. With 4% of the total SN en- 
ergy in the form of kinetic energy input as suggested by 
Dorfi (1993) would raise the total star formation rate to 
M* ~ 3.3 M yr _1 . 

7. Application to the Galaxy 

In this Section we will investigate how close the simple 
analytical models of Sect. 14.21 and IB1 can get to the obser- 
vations of the Galaxy. 

7.1. Self-gravitating disk in z direction (SGZ) 

We use Q = 1.0, S = 5, a = 3 10 7 yrM pc- 3 , M = 
10 -1 M yr" 1 , and £ = 4.6 1CT 8 (pc/yr) 2 . Fig.Hshows the 
derived turbulent velocity, driving length scale, gas surface 
density, disk scale height, gas density, viscosity, molecular 
fraction, volume filling factor, and star formation rate for 
three different rotation velocities. For the model of the 
Galaxy we adopt v Iot — 220 kms -1 . This results in a 
total gas mass of M gas ~ 6.8 10 9 M , a total molecular 
mass of M mo \ ~ 3.3 10 9 M , a total atomic gas mass of 
Mhi ~ 3.5 10 9 M Q , and a total star formation rate of 
M* ~ 2.6 M yr _1 . 

The local density p ~ 2 cm -3 and the local surface 
density of £ ~ 2 10 21 cm -2 at the solar radius (Binney & 
Trcmainc 1987) are well fitted by our model. At the solar 
radius the disk scale is H ~ 400 pc, the turbulent driving 
scale length is ^driv ~ 90 pc, and the effective viscosity is 
v ~ 27 10 25 cm 2 s _1 . The driving wavelength is compara- 
ble to the generally accepted value of 50 pc< /driv <150 pc 
(see e.g. Ruzmaikin et al. 1988). Moreover, the derived star 
formation rate and gas surface density are comparable to 
those observed for the Galaxy. 

The constant volume filling factor of <j>y ~ 2 10~ 3 
gives an equivalent Reynolds number of the fully grav- 
itational model (Vollmer & Beckert 2002) of Re = 22 
for a non-Kolmogorov spectrum (D = 2). In Vollmer & 
Beckert (2002) we used Q = 1, M = 10~ 2 Iv^yr" 1 , 
a = 10 7 yrM pc -3 and Re = 50. All these parameters 
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Fig. 4. Radial profiles of disk parameters for the case of 
the SGZ model for three different rotation velocities (solid 
lines: u ro t=250 kms -1 , dashed lines: 175 kms' 1 , dotted 
lines: 100 kms -1 ). 



are within a factor 2 consistent with the values described 
above. This again shows the equivalence of the two mod- 
els. 

7.2. SN remnant size as driving length scale (SNR) 

For this model we use Q = 1, S = 5, a = 2 10 7 yr M Q pc~ 3 , 
M = 10- 1 Moyr" 1 , and £ = 4.6 10 -8 (pc/yr) 2 . The de- 
rived disk properties can be seen in Fig. For a rotation 
velocity of u ro t = 220 kms -1 , we obtain the following re- 
sults: total gas mass M gas ~ 6.3 10 9 M Q , total molecular 
mass of M mo \ ~ 2.0 10 9 Mq, a total atomic gas mass of 
Mhi ~ 4.3 10 9 M , and the total star formation rate given 
in Sect. El 

The driving wavelength at the solar radius is 
i(jriv=100 pc. The other local quantities (S, p, H) are 
comparable to those of Sect. 17.11 As already mentioned 
the star formation rate calculated using Eq. (|20|> is some- 
what too high and that calculated using the energy flux 
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Fig. 5. Radial profiles of disk parameters for the case of 
the SNR model for three different rotation velocities (solid 
lines: t> ro t=250 kms -1 , dashed lines: 175 kms -1 , dotted 
lines: 100 kins" 1 ). 

conservation (Eq. is somewhat too low compared to 

the observed value. 

8. A realistic gravitational potential 

In this Section we use an analytic gravitational potential 
$ g given by Allen & Santillan (1991). It consists of a disk, 
bulge, and halo component. We use the following values 
for the parameters of Allen & Santillan (1991): 

- bulge: Mi = 1.41 10 10 M Q , b x = 387 pc, 

- disk: M 2 = 8.56 10 10 M Q , a 2 = 5318 pc, b 2 = 250 pc, 

- halo: Mi = 1.07 10 11 M Q , a 3 = 12 kpc. 



P«turb 



Vertical pressure equilibrium is described by 

dz 



1 -S(tGE + ^ 



(76) 



8.1. Self-gravitating disk in z direction 

We adopt Q = 1.0, 5 = 5, M = 5 10~ 2 Meyr" 1 , and 
a = 410 7 yrM pc -3 . 




Fig. 6. SGZ model. Solid lines: radial profiles of disk 
parameters for a realistic gravitational potential. Dashed 
lines: disk properties including magnetic and thermal pres- 
sure. 

The resulting radial profiles for the disk parameters 
are shown in F'ig. EI It turns out that the gas disk is close 
to the stage of self-gravitation in z direction, but is still 
dominated by the gravitational potential of disk, bulge, 
and halo. The total gas mass is M to t = 7.5 10 9 M©, the 
total molecular mass is M mo i = 2.8 10 9 M©, and the star 
formation rate is Af» = 5.0 M© yr _1 . 

For R > 5 kpc the gas surface density has approxi- 
mately a 1/ R profile. The disk is flaring with an increasing 
flaring angle for increasing R. The volume filling factor is 
still approximately constant with radius. The star forma- 
tion rate can be described by an exponential with a scale 
length of ~3.5 kpc for radii between 3 and 10 kpc. 

8.2. SN remnant size as driving length scale (SNR) 

In the case of the SNR model we adopt slightly dif- 
ferent parameters in order to match Galactic observa- 
tions: Q = 1.0, <5 = 5, M = 510~ 2 Moyr" 1 , and 
a = 2 10 7 yrM©pc~ 3 . The resulting radial profiles for 
the disk parameters are shown in Fig. The total gas 
mass is M to t = 6.5 10 9 M©, the total molecular mass 
is M mo i = 2.7 10 9 M©, and the star formation rate is 
M* = 6.5 M©yr _1 . All disk properties have values com- 
parable to those of the SGZ model. The main differences 
to the SGZ model are that (i) the volume filling factor 
increases significantly with decreasing galactic radius, (ii) 
the star formation rate rises more steeply to the galaxy 
center, and (iii) the molecular fraction rises less steeply 
to the galaxy center than for the SGZ model. As for the 
SGZ model the gas surface density has approximately a 
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Fig. 7. SNR model. Solid lines: radial profiles of disk 
parameters for a realistic gravitational potential. Dashed 
lines: disk properties including magnetic and thermal pres- 
sure. 




1 o i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

2 4 6 8 10 12 14 

radius (kpc) 

Fig. 8. Radial profiles of the total (dotted line), molecular 
(solid line), and atomic (dashed line) gas surface density. 

l/R profile for R > 5 kpc. The star formation rate shows 
approximately a l/R 2 profile. 

8.3. The Hi to CO surface density relation 

Fig. [S] shows the radial profiles of the total, molecular, 
and atomic gas surface density fro the SGZ model. These 
profiles well resemble those found observationally by Wong 
& Blitz (2002). The authors give the following relation 
between the atomic and the molecular gas surface density: 
Shi/Sco oc R 15 . In order to check our model for this 




1 .0 r 
0.5 r 
0.0 L 
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radius (kpc) 

Fig. 9. The radial dependence of the ratio between the 
atomic and the molecular gas surface density Ehi/Sco 
(upper graph) and its slope (lower graph). 

relation we show in Fig. the radial dependence of the 
ratio Sri/Sco and its slope for the model SGZ. The curve 
for the SNR model has approximately the same slope. For 
R > 3 kpc both models reproduce the observed exponent 
of 1.5. 

8.4. Magnetic fields and thermal pressure 

For the vertical pressure equilibrium (Eq. i|13|) . HMD ) 
we have neglected the magnetic and thermal pressure: 
Pmagn = B 2 / (8ir) and pth — n fee T, where B is the mag- 
netic field strength, n is the particle density, fee is the 
Boltzmann constant and T is the gas temperature. 

We have solved the set of equations of Sect. |S] including 
these pressure components and setting the turbulent pres- 
sure to pturb = \ P ^t ur b- We assume equipartition between 
the energy density of the turbulent gas and the magnetic 
field \pv 2 mh = B 2 /(8n) and T^IOOOO K. The resulting 
disk properties are plotted as dashed lines in Fig.[f)]and[7| 

The magnetic field derived from equipartition de- 
creases from ~20 pG in the center to ~1 pG at the edge 
of the disk. Its value at 8.5 kpc is ~7 pG. Due to the 
additional pressure terms the turbulent velocity, gas col- 
umn density, gas density, molecular fraction, and the star 
formation rate decrease, whereas the driving length scale, 
disk scale height, and the viscosity increase. The volume 
filling factor decreases slightly in the case of the SGZ 
model and increases in the case of the SNR model. 

The star formation rates are M» = 3.6 M Q yr _1 for the 
SGZ model and M* = 4.2 Moyr" 1 for the SNR model. 
The radial profile of the magnetic field then fits nicely that 



12 



B. Vollmer and T. Beckert: Turbulent viscosity in clumpy accretion disks II 




10 

gas surface density (M Q pc~ 2 ) 



Fig. 10. The star formation rate E» as a function of the 
total gas surface density E. Solid line: SGZ model. Dashed 
line: SNR model Dotted lines: E* oc E 1 - 5 and E* oc E 2 . 

proposed by Strong, Moskalenko, & Reimer (2000) in its 
absolute values. It can be described by an exponential with 
a scale length of 8 kpc and a local value of the magnetic 
field of ~7 /zG for both models. 

8.5. The star formation law 

For the relation between the star formation and the total 
gas density, Wong & Blitz (2002) used a Schmidt law of 
the form 

E* cx E^ . (77) 

They found /3 = 1.7 ± 0.3 if they use the star formation 
derived from Ha measurements corrected by an extinction 
that is proportional to E. Fig.llOlshows the star formation 
rate E* as a function of the total gas surface density E 
for our model together with two curves for = 1.5 and 
(3 = 2. The star formation rate of the SGZ model has 
an exponent (3 = 2 for E < 50 M Q pc -2 and /3 = 1.5 
for E > 50 M Q pc~ 2 . The star formation rate of the SNR 
model has always an exponent greater than 2 and steepens 
considerably for E > 50 M Q pc~ 2 . Thus, only the SGZ 
model can reproduce the observational results of Wong & 
Blitz (2002). 

9. Discussion 

The dominating stellar disk (DSD) model applies for gas 
deficient galaxies (SO and cluster galaxies). These are sep- 
arate cases and will be investigated in more detail in a 
separate paper. 

The model of a self-gravitating gas disk (SGZ model) 
and the model where the SN remnant size is assumed to 
be the driving length (SNR model) can be successfully 
applied to the Galaxy. The input parameters are Q ~1, 
(5 = 5, M ~ 0.05 - 0.1 Moyr- 1 , £ ~ 4.6 10~ 8 (pc/yr) 2 , 
and a ~ 2 — 3 10 7 yrM©pc~ 3 . The constants of molecule 



formation a in all models are very close to that given by 
Hollenbach & Tielens (1997). The fully analytical models 
(Sect. 14.21 and EJ) both describe the gas properties and 
the star formation rate of the Galaxy well. Their radial 
dependencies are very similar. 

The models with a realistic gravitational potential fit 
observations well and lead to star formation rates that 
are consistent with the observed value. The star forma- 
tion rate in the inner part of the galactic disk is more 
realistic for the SGZ than for the SNR model. Moreover, 
only the SGZ model reproduces the observed exponen- 
tial radial profile of the star formation rate with the right 
length scale. Thus, the fit of the radial profile of the star 
formation rate favors the SGZ model. 

Including a magnetic field in equipartition with the 
turbulent energy density leads to a realistic radial profile 
of the magnetic field in accordance with derived profiles 
in the literature. 

The main difference between the two models is found 
in the radial profile of the volume filling factor of self- 
gravitating clouds, which is not well constrained by obser- 
vations. Blitz (1993) states that a "glaring deficiency in 
galactic studies of giant molecular clouds is a quantitative 
study" of the radial dependence of their properties. Thus, 
both models can be regarded as equally valid. 

We conclude that our models correctly reproduce the 
following observations: 

— the integrated total, molecular, and atomic gas mass, 

— the integrated star formation rate within a factor of 2, 

— the radial dependence of the ratio between the atomic 
and molecular gas surface densities. 

The exponent for the Schmidt law might be ~20% too 
high. 

In Sect. we show that the viscosity and star forma- 
tion prescriptions for the SN driven and the gravity-driven 
turbulence (without SN feedback) are equivalent for Q = 1 
and 6=1. However, the absolute values of the energy flux 
-£§j = -Ei/^n- (Eq. (ESI) is different for the two mod- 
els. In the case of gravity-driven turbulence the energy flux 
conservation equation together with the angular momen- 
tum equation and the vertical pressure equilibrium lead to 
Zdriv = H independent of a specific viscosity prescription. 
For the viscosity v = 7Uturb^driv (Eq. ©) the gravity- 
driven turbulence model reads 7 = Re^ 1 (Eq. H59|) ~) and 
H/ldiiv — 1, whereas the SGZ model gives 7=1 and 
H/ldxiv — 4-5. Thus the energy flux due to SN is a factor 
of ^200 larger than that due to gravitational instabilities. 
This implies that whenever the star formation exceeds 1% 
of its equilibrium value for a given galactic gas disk (see 
Sect. l4.2|) . the energy input due to SN dominates over that 
due to gravitational instabilities. Whenever it is smaller 
than 1% of its equilibrium value, gravitational instabili- 
ties take over and drive the ISM turbulence. In this case 
the mass accretion rate M is ~10 times smaller than that 
of the SN driven turbulent gas disk. 
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10. Star formation 

In this section we study the behaviour of the star for- 
mation efficiency (SFE) defined as the inverse of the star 
formation time scale: 

SFE = (U)- 1 = ^ . (78) 

We are especially interested in the dependence of the SFE 
on Q. 



dependent equation of transport of angular momentum. 
Thus, this estimate must be regarded with caution. Our 
equilibrium model suggests that at the beginning of a star- 
burst the turbulent velocity stays constant at 10 kms^ 1 , 
whereas the driving length scale is reduced. This coun- 
terbalances the enhanced star formation rate. Evidently, 
after a short time the available gas will be consumed and 
the SNe might heat the disk and lead rapidly to a Q > 1 
disk. 



10.1. Starbursts (Q < 1) 

There are two possible definitions of a starburst: 

1. definition based on the absolute value of the star for- 
mation rate measured by the Ha, radio continuum, 
or FIR. In this case, the star formation is enhanced, 
because of an enhanced gas surface density, e.g. the 
Schmidt law (£* oc E' 3 ) is still valid. 

2. definition based on the SFE. In this case the SFE is 
enhanced (the star formation time scale is reduced). 

Rownd & Young (1999) studied a sample of starburst 
galaxies. Out of 85 galaxies only 24 are starburst galaxies 
according to definition 2. Wong & Blitz (2000) observed 
one of them, NGC 4736. They found, in agreement with 
Rownd & Young (1999), that the nuclear Ha ring does 
not follow a Schmidt law, but has a much higher SFE 
than expected by a Schmidt law. 
Our SGZ model yields 

SFE = cx Q- 1 (79) 

Zj 

for the SGZ model and 

SFE = cx Q- 177 (80) 

for the SNR model. In Sect. 14.21 we have shown that for 
the SGZ model the turbulent velocity dispersion is inde- 
pendent of Q (Eq. lHHQ an d thus Q measures the ratio 
between the total enclosed mass and the total gas mass, 
so that for a given galactic potential Q cx E _1 . This is 
also approximately valid for the SNR model. 

If the infall time scale is short enough that a consid- 
erable fraction of gas can be transported to small galactic 
radii where it is accumulated, star formation will be en- 
hanced much stronger than predicted by the Schmidt law 
in this region. This is one possibility to create a central 
starburst. The massive infall needed might be due to a 
tidal interaction with another galaxy, an interaction be- 
tween the ISM and the intracluster medium in a galaxy 
cluster, or a bar instability. Based on Eq. and (|S7Jl 
our model would then predict 

£* cx S 3 . (81) 

However, one has to keep in mind that a Q < 1 disk is 
not a stable configuration and can in principle not be de- 
scribed by an equilibrium disk. One has to solve the time 



10.2. Exhausted star formation (Q > I) 

In the opposite case, when gas is taken away from the 
galaxy by enhanced star formation or external effects, Q 
becomes greater than unity. Then the DSD model applies 
giving 

SFE = — ^ocQ~s (82) 

for the SGZ model. In this case star formation becomes 
less efficient, but this effect is not as pronounced as the 
increase for Q < 1 discussed above in Sect. 110.11 

11. Summary and Conclusions 

We extend the model of a turbulent cloudy gas disk 
(Vollmer & Beckert 2002) by including energy input 
through SN explosions. This is realized (i) by adapting 
the energy flux conservation equation where the energy 
flux transferred by turbulence to smaller length scales is 
balanced by the energy input due to SN explosions (SGZ 
model) or alternatively (ii) by assuming that the size of 
a SN remnant equals the turbulent driving length scale 
(SNR model). In the SGZ model the SN energy flux is as- 
sumed to be proportional to the local star formation rate. 
The local star formation rate p» is assumed to be propor- 
tional to the mean density and inversely proportional to 
the local free fall time of the clouds. The factor of propor- 
tionality is the probability to find a self-gravitating cloud, 
i.e. the volume filling factor. The integration length in z di- 
rection is assumed to be the turbulent driving length scale, 
i.e. the length scale over which clouds are self-gravitating: 

£* P* ^driv- 

We assume the turbulence to be intermittent due to 
self-gravity. The turbulence v is thus reduced to v = 
7^turb ^driv with 7 < 1 . For Q < 1 we find 7=1. 

For model (i) we have calculated two separate cases: 

— a dominating stellar disk mass (DSD) (Q > 1) and 

— a self-gravitating gas disk in z direction (SGZ) (Q ~ 
!)• 

For the SGZ model and tt' ~ —ft/R, the radial depen- 
dences of the disk properties (height, turbulent driving 
length scale, surface density, density, viscosity, volume fill- 
ing factor, star formation rate, and molecular fraction) can 
be expressed as functions of the mass accretion rate M, 
Q, the cloud size divided by the driving length scale 5, 
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the fraction of SN energy that is transformed into kinetic 
energy £, the constant of molecule formation a, and the 
rotation velocity £l(R). In the case of the DSD model, the 
dependences on f2(i?) are the same as for the SGZ model, 
whereas those on M, Q, a, and £ are different. In addi- 
tion, a factor that depends on the rotation velocity w rot 
appears. 

The SNR model gives radial dependences of the disk 
properties with exponents that are close to those of the 
SGZ model. Both analytical models (SGZ and SNR) re- 
produce an almost constant turbulent gas velocity as ob- 
served. 

For the SGZ model we conclude that 

1. the ratio between the turbulent driving length scale 
^driv and the disk height H equals the square root of 
the volume filling factor (f>y, 

2. for Q = 1 the viscosity and star formation prescrip- 
tions are equivalent to those given for the fully gravi- 
tational model (Vollmer & Beckert 2002), 

3. the radial dependences are the same as those for the 
fully gravitational model. 

The SGZ and SNR models are applied to the Galaxy 
and give a good description of its gas disk. We derive Q <~ 
1 and a mass accretion rate within the disk of M <~ 0.05 — 
0.1 Moyr" 1 , which is a factor 5-10 higher than that of 
the fully gravitational model. 

We include a realistic gravitational potential into the 
model, which improves the fit to observations. In this case 
the SGZ and SNR models give equally good fits to obser- 
vations. 

The influence of magnetic fields and thermal pressure 
is investigated. If both are included in the model, the tur- 
bulent velocity decreases with respect to the case where 
only the turbulent pressure is taken into account. An in- 
crease of the mass accretion rate by a factor 1.5 makes the 
turbulent velocity dispersion increase again to 10 kms -1 . 
In this case the radial profile of the magnetic field is close 
to observations. 

The models nicely reproduce the observed radial de- 
pendence of the ratio between the Hi and molecular sur- 
face densities. We derive a Schmidt law of the form £ cx 
E 2 . A possible explanation for the enhanced star forma- 
tion efficiency observed in a few spiral galaxies is given. 

We conclude that for Q = 1 the models for turbu- 
lence (i) driven by SN energy input and (ii) generated by 
instabilities involving self-gravitation and maintained by 
energy input from differential rotation and mass transfer 
are consistent. The main difference in the disk properties 
is that Idiw/H < 1 for case (i) and l^v/H = 1 for case 
(ii). 

Both mechanisms are not exclusive and might possibly 
coexist in galactic disks. A galaxy that slowly forms stars 
will consume its gas and run into the regime Q > 1. Then, 
the timescale for star formation will increase. If locally 
star formation stops entirely, the turbulent energy will be 
dissipated within /driv/wturb ~ 10 7 yr, the disk shrinks in 
z direction and ultimately becomes self-gravitating in z. 



This is the point, where the fully gravitational model of 
Vollmer & Beckert (2002) applies. In this case the driving 
length scale increases and equals the disk height. Thus, 
the energy dissipation rate, which is proportional to l^? v 
decreases and the small energy supply due to differential 
rotation is large enough to maintain turbulence. 

Such a disk has Q > 1 and a 5 times lower mass ac- 
cretion rate than a star forming disk. Once the gas disk is 
self-gravitating in z it will form stars again and switch to 
the SGZ/SNR state. 

The occurrence of a bar or an external accretion event 
of a companion can lead to a temporal increase of the 
mass accretion rate. If Q becomes smaller than 1 during 
this event, the star formation timescale decreases rapidly. 
The galaxy will then again reach the regime of Q ~ 1 . The 
detailed evolution of such a process can only be investi- 
gated by solving the time dependent equation for angular 
momentum transport. 
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